clear

c1 = rand(3,1);
c2 = rand(3,1);
C1 = expm(so3(c1));
C2 = expm(so3(c2));

theta12 = vee(logm(C1'*C2))
theta21 = vee(logm(C2*C1'))